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background: Cytosine arabinoside-based chemotherapy coupled with anthracycline is currently the first-line treatment for acute 
myeloid leukaemia (AMI), but diverse responses to the regimen constitute obstacles to successful treatment. Therefore, outcome 
prediction to chemotherapy at diagnosis is believed to be a critical consideration. 

methods: The mRNA expression of 1 2 genes closely involved in the actions of cytosine arabinoside and anthracycline was evaluated 
by real-time reverse transcriptase PGR (RT-PCR), in 54 diagnostic bone marrow specimens of M2-subtype AMI. 
results: Low expression levels of ribonucleotide reductase MI (RRMI) and high expression levels of topoisomerase 2 beta (T0P2B) 
were correlated with longer survival in a univariate analysis. Another interesting finding is that high ratios of TOP2B/RR/V\2 and T0P2BI 
T0P2 alpha (T0P2A) in a combined analysis were also shown to have a prognostic impact for longer sun/ival with improved accuracy. 
Among the four markers, when adjusted for the influence of other clinical factors in multivariate analysis, the TOP2B/TOP2A ratio was 
significantly correlated with treatment outcomes; patients with high ratios trended toward longer disease-free sun/ival (HR, 0.24; 
P = 0.002) and overall survival (HR, 0.29; P = 0.005). 

conclusion: Genes with distinct expression profiles such as TOP2B/TOP2A expression ratio at diagnosis can be employed for 
outcome prediction after the treatment with standard regimens in AML patients with M2 subtype. 
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Acute myeloid leukaemia (AML) is a heterogeneous haematologic 
malignancy and is the most common type of the condition, 
accounting for up to 60% of leukaemia cases. The French- 
American-British classification subdivides AML into nine sub- 
types (M0-M7 and M4Eo), which can be distinguished by their 
lineages and the degree of differentiation of myeloid cells 
(Lowenberg et al, 1999). Among these, the M2 subtype is the 
top-ranked disease ( ~ 30%) in terms of incidence. 

Nucleoside analogues or combination with anthracyclines have 
been used for the majority of AML patients and are generally 
regarded as the most effective therapies (Wiernik et al, 1992). 
Although approximately two-thirds of the patients, at least, achieve 
complete remission after receiving intensive induction regimens, 
most of these patients suffer from relapsed diseases. Thus, a failure 
of first-line chemotherapy in these patients can also attenuate the 
antitumour effect of subsequent chemotherapeutics, owing to 
substantial multi-drug resistance acquisition (Greenberg et al, 
2004). At diagnosis, therefore, it is believed that the accurate 
response prediction against initial regimens is an extremely 
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important consideration, and can provide another option for 
AML treatment. 

Once cytosine arabinoside, a pyrimidine analogue, enters into 
the cells through specialised transporters such as human 
equilibrative nucleoside transporters (hENT) and concentrative 
nucleoside transporters (hCNT), it is finally converted into an 
active drug by deoxycytidine kinase (dCK) (Stegmann et a/, 1995; 
Hubeek et al, 2005). Recent studies using clinical samples have 
demonstrated that high levels of dCK and hCNT3 expression were 
correlated with favourable responses among patients with AML 
and CLL, respectively, who received nucleoside drug-based 
chemotherapy (Mackey et al, 2005; Yamauchi et al, 2009). 
Conversely, a reversal of drug metabolism by ecto- 5' -nucleotidase 
(CD73), 5',3'-nucleotidase (NT5C) and cytidine deaminase (CDA) 
has been regarded as indicators of poor outcome in AML 
treatment (Ujhazy et a/, 1996; Galmarini et al, 2004; Mahlknecht 
et aly 2009). Anthracyclines such as idarubicin and daunorubicin, 
supplemental drugs used in the treatment of AML, exert cytotoxic 
effects by intercalating into DNA strand and inhibiting topoi- 
somerase activities in cancer cells (Carella et al, 1990; Richardson 
and Johnson, 1997). Although there have been attempts to 
elucidate the relationship between topoisomerase II expression 
and drug response in clinical AML samples (Kaufmann et al, 1994; 
McKenna et al, 1994; Wang et al, 2009), it still remains unsettled. 



Another drug resistance mechanism in the chemotherapy contain- 
ing nucleotide analogues is the reduction of intracellular drug 
accumulation by ATP-binding cassette (ABC) transporters (Pogliani 
et aly 2001; Steinbach et a/, 2006). These findings imply that the genes 
encoding for molecules associated with action mechanism and 
delivery of the induction regimens are potential predictive parameters 
of clinical outcomes in the standard chemotherapy for AML. 

In this study, in an effort to address the clinical significance of 
12 candidate genes that perform critical functions in the actions of 
cytosine arabinoside and anthracycline, we evaluated their mRNA 
levels in bone marrow (BM) samples from 54 diagnostic AML 
patients with the M2 subtype. The findings of correlation analyses 
of gene profiles with clinical outcomes may provide us with 
valuable insights into possible strategies to improve treatment 
outcomes in AML therapy. 



MATERIALS AND METHODS 
Clinical samples 

BM aspirates of 54 adults newly diagnosed with the M2 subtype of 
AML and 4 healthy donors were collected from Catholic University 
Seoul St. Mary's Hospital and Chonnam National University 
Hwasun Hospital (Table 1). Written informed consents for the 
cryopreservation and use for further research of all BM samples 
were obtained from the patients. The Institutional Review Board of 
Catholic University Seoul St Mary's Hospital and Chonnam 
National University Hwasun Hospital approved all research on 
the human subjects participating in this study. The patient 
population consisted of 28 males and 26 females with a median 
age of 42.5 years (range, 20-72). Induction chemotherapies were 
executed by using a combination of l-j5-D-arabinofuranosyl 
cytosine (Ara-C, Ar= 14) or N4-behenoyl-l-j5-D-arabinofuranosyl 
cytosine (BH-AC, N= 35) combined with idarubicin. Four patients 
received Ara-C alone and one did Ara-C with etoposide. Among 
these, 33 patients (61.1%) received haematopoietic stem cell 



Table 



Patient characteristics 



Category 



Value 



Total specimens, no. 

Age, median year (range) 
Sex, no. (%) 

Male 

Fennale 

BM blasts, median % (range) 

FAB subtype, no. (%) 
M2 

Others 

Cytogenetics, no. (%) 
Normal 
Abnormal 

t(8;2l) 

Others^ 

HSCl no. {%) 
Yes 
No 

Disease-free survival, day (range) 
Overall survival, day (range) 



54 

40 (20-72) 

28 (5 1 .8) 
26 (48.2) 

42.5 (2-92) 



54 (100) 
0(0) 



1 6 (29.6) 

38 (70.4) 

24 (44.4) 

14 (26.0) 



33 (61.1) 
21 (38.9) 

0-1370 
9-1510 



Abbreviations: BM = bone marrow; HSCT = haematopoietic stem cell transplanta- 
tion. "Other cytogenetics include t( I ;4), + 8/t( I ; I I )/t(6;9)/t( I I ; 1 2)/t( 1 2;22)/del(5q)/ 
del(9q)/del( 1 6q)/add(2)/inv( 16)/ - 1 0, - 1 3, - 1 5, - 1 9/ + 8/ + 1 4, + 20/ 1 3 
cenh-K abnormalities. 
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transplantation (HSCT) after induction and consolidation chemo- 
therapies. A total of 16 patients (29.6%) evidenced no cytogenetic 
changes, and the others (70.4%) exhibited cytogenetical abnorm- 
alities including t(8;21). Additionally, normal BM samples from 
four healthy donors were used as references throughout the 
experiments and data analyses. 

RNA preparation and real-time PGR 

All of the BM mononuclear cells (MNCs) from patients and from 
healthy donors were kept at — 80 °C until use. Total RNA extracted 
using a Qiagen RNA isolation kit (Hilden, Germany) was converted 
to cDNA with a reverse transcription kit (Invitrogen, Carlsbad, CA, 
USA). The cDNA was mixed with SYBR Green PGR Master Mix (PE 
Applied Biosystems, Foster Gity, GA, USA) and the specific 
primers (Supplemental Table SI) for 12 genes including CDl^^ 
NT5C, CD A, dCK, hCNT3, hENTl, ribonucleotide reductase Ml 
(RRMl), RRM2, topoisomerase 2 alpha {T0P2A), T0P2 beta 
(T0P2Bi ABC subfamily A member 3 (ABCA3) and ABCBL The 
candidate genes and GAPDH were amplified in triplicate using an 
ABI Prism 7900 Sequence Detection System (PE Applied 
Biosystems) under the following conditions: 40 cycles of 95 °G 
for 30 s, 60 °G for 30 s and 72 °G for 30 s. To normalise the raw data, 
ACt was generated by subtracting the Ct value for GAPDH from 
the Ct of 12 genes in an individual patient sample. The results were 
subsequently adjusted to the average gene levels in four healthy 
donors (control ACt — patient ACt = AACt). The results were 
expressed in arbitrary quantitative PGR units (2~^^^^^^). In some 
cases, these data were applied to gene-clustering analysis to 
explore whether there was any unique pattern between the 
individual gene expression. 

Correlation analysis of gene expression with 
clinical outcomes 

The mean expression values for each gene in four healthy donors 
were set as the cutoff point 1.0 and the gene levels in patient 
samples were expressed as relative values compared with the mean 
values of each gene from healthy donors. We then dichotomised 
total study population of 54 members into two groups for 12 
individual genes by the cutoff value. Populations with expression 
level lower and higher than 1.0 were named the low and high 
groups, respectively. We defined DPS and OS as the time from first 
remission to relapse or death and the time of diagnosis to death or 
date of last follow-up, respectively. Using both populations clustered 
according to the levels of each gene and their clinical information, 
Kaplan-Meier estimate constructed survival curves for DPS and OS. 
Differences in the clinical outcomes and statistical significance 
between both groups were then determined using a univariate 
analysis with log-rank test. In combined analysis using the expres- 
sion ratios of the two genes, especially, the patients were grouped 
by the median value of the gene ratio. Cox proportional hazards 
model was employed to investigate the influence of clinical factors 
against the prognostic gene markers. The details of this are provided 
in the figure legends. All statistic analyses were done by MS Office 
Excel program and SPSS 12.0 software (Chicago, IL, USA). A level of 
P<0.05 was regarded as statistically significant in this study. 

Cell lines and maintenance 

Human leukaemia cell OCI/AML-2 (AML-2) was cultured in alpha- 
MEM media (Gibco, Grand Island, NY, USA) supplemented with 
10% PBS and antibiotics. We generated and maintained AML-2/ 
ID AC cells double-resistant to Ara-C and idarubicin as previously 
described (Song et al, 2009a). 

Methylation- specific PGR (MSP) analysis 

Cell lines were rinsed twice in PBS, and their gDNA were isolated using 
AccuPrep Genomic DNA Extraction kit (Bioneer, Daejeon, Korea). 
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Sodium bisulfite- converted DNA was then obtained from 0.5 fig of 
the gDNA using an EZ DNA Methylation kit (Zymo Research, 
Irvine, CA, USA). Total DNA of 15 ng were used for MSP 
corresponsive to CpG islands. The primer pairs were designed 
using MethPrimer software (Li and Dahiya, 2002), and are as 
follows: RRM2, Fw: 5'-ttagtttggtcgataaggagaatac-3' and Re: 5'-taa 
aaaaacaatacgactttaacgc-3' for methylation (305 to 566 from 
transcription start site, TSS), and Fw: 5'-tagtttggttgataaggagaata 
tgg-3' and Re: 5'-aaaaaaacaatacaactttaacacc-3' for unmethylation 
(306 to 565 from TSS); T0P2B, Fw: 5'-gatagtatggcgttgattgatagc-3' 
and Re: 5'-tacttcaaaaacaaccttaacctcg-3' for methylation ( — 679 to 
— 476 from TSS), and Fw: 5'-tggatagtatggtgttgattgatagt-3' and Re: 
5'-cttcaaaaacaaccttaacctcact-3' for unmethylation ( — 681 to —78 
from TSS); T0P2A, Fw: 5'-ggttggagtgtagtggagtaattac-3' and Re: 
5'-aacacattaaaaaaccgaaacga-3' for methylation ( — 792 to — 64 from 
TSS), and Fw: 5'-aggttggagtgtagtggagtaattat-3' and Re: 5'-caacacat 
taaaaaaccaaaacaaa-3' for unmethylation ( — 793 to — 63 from TSS). 
The products were analysed by agarose gel electrophoresis. 

Semi-quantitative RT-PCR analysis 

Total RNA from cell lines was converted to cDNA using a 
CycleScript Reverse Transcriptase kit (Bioneer). After mixing 
the cDNA with specific primers for RRM2, T0P2A and T0P2B, 
PGR was carried out. GAPDH was used as an internal control. 
The product band intensity was estimated using Image J software 
(http://rsbweb.nih.gov/ij/) and normalised using GAPDH. The 
quantitative results were then displayed in terms of relative 
expression to the wild-type control cells. The result analysis was 
performed using Student's t-test. 

RESULTS 

Patient population and clinical features 

We obtained 54 BM MNCs from patients with the M2 subtype of 
AML at diagnosis; their details are described in Table 1. A total 
of 16 patients exhibited normal karyotype and 24 patients had 
t(8;21). Each of the others evidenced cytogenetically distinct 



abnormalities. For remission induction, a majority of the patients 
received standard chemotherapy with Ara-C (Ar= 14) or BH-AC 
(Ar=35) coupled with idarubicin. Four of the remaining patients 
were treated with Ara-C alone, and one received a regimen of Ara- 
C plus etoposide. After induction and subsequent consolidation 
chemotherapies, 33 patients (61.1%) received HSCT. Median DFS 
and OS periods for all populations were 521 (95% CI, 99-943) and 
737 (95% CI, 222-1252) days, respectively. In a univariate analysis, 
we noted that HSCT resulted in prolonged patient survival 
(P<0.05). Patients with t(8;21) also evidenced better responses 
than the population with normal or other cytogenetics. The 
correlation between other clinical features (age, sex, RI regimen 
and blast % in BM) and treatment outcomes did not reach 
statistical significance at the specified level (Table 2). 

Transcriptional levels in AML samples and their 
correlation with clinical outcomes 

Expression profiles for 12 genes were analysed in a total of 54 
samples to identify a prediction marker for treatment outcome at 
diagnosis. Relative expression levels from real-time PCR analysis 
of 12 genes were displayed in box plots (Figure IB). In order to 
characterise the relationship of the gene expression pattern with 
clinical outcomes, we clustered the patients into two groups 
according to the mRNA levels of each of the genes relative to the 
levels of healthy control subjects. The prognostic impacts of 
individual genes between the low and high groups were evaluated 
via a univariate analysis by Kaplan-Meier estimation and log- rank 
test, and summarised (Table 2). In this analysis, we noted that the 
high RRM2-level group (N=18) evidenced 2.85-fold (95% CI, 
1.22-6.62; P = 0.002) and 2.69-fold (95% CI, 1.17-6.20; P = 0.005) 
increases in the hazard ratio (HR) for DFS and OS, respectively, 
relative to the low group (Figure 2). High T0P2B at the 
transcriptional level was also associated with longer DFS 
(P = 0.023) and OS (P = 0.024) in the AML patients. Because the 
two clinical factors, cytogenetics abnormality and SCBT, already 
evidenced an association with treatment outcome, we investigated 
whether there was a correlation between both factors and 
biological markers (Supplemental Table S2). In the study, no 



Table 2 Univariate analysis of treatment outcome using 1 2 


genes and clinical factors by log 


;-rank test 






Variables 


Group (N) 


Disease-free survival 


Overall survival 




HR (95% CI) 


P-value 


HR (95% CI) 


P-value 


CD73 


H (21) vs L (33) 


0.77 (0.36-1.64) 


0.500 


0.71 (0.34-1.50) 


0.391 


NT5C 


H (48) vs L (6) 


0.57 (0. 1 3-2.59) 


0.322 


0.67 (0.20-2.23) 


0.509 


CDA 


H (6) vs L (48) 


1.29 (0.40-4.1 1) 


0.630 


1 . 1 1 (0.37-3.33) 


0.842 


dCK 


H (33) vs L (21) 


1.07 (0.51-2.26) 


0.965 


0.98 (0.46-2.09) 


0.852 


hCNT3 


H (39) vs L (15) 


0.73 (0.31-1.70) 


0.416 


0.69 (0.29-1.63) 


0.349 


hENTI 


H (36) vs L (18) 


1 .75 (0.80-3.83) 


0.202 


1.25 (0.53-2.90) 


0.621 


RRMI 


H (1 1) vs L (43) 


1.62 (0.61-4.32) 


0.247 


1.56 (0.59-4.10) 


0.303 


RRM2 


H (18) vs L (36) 


2.85 (1.22-6.62) 


0.002 


2.69 (1.17-6.20) 


0.005 


TOP2A 


H (26) vs L (28) 


1.48 (0.71-3.07) 


0.278 


1 .29 (0.62-2.68) 


0.487 


TOP2B 


H (15) vs L (39) 


0.36 (0.17-0.76) 


0.023 


0.35 (0. 1 6-0.74) 


0.024 


ABCA3 


H (52) vs L (2) 


0.41 (0.05-3.53) 


0.187 


0.57 (0.09-3.61) 


0.430 


ABCBI 


H (30) vs L (24) 


0.78 (0.37-1.62) 


0.484 


0.73 (0.35-1.53) 


0.386 


Age 


^60 (13) vs <60 (41) 


1 .87 (0.75-4.70) 


0.099 


1 .94 (0.77-4.92) 


0.090 


Sex 


Female (26) vs Male (28) 


0.63 (0.30-1.32) 


0.201 


0.69 (0.33-1.43) 


0.316 


HSCT 


No (21) vs Yes (33) 


2.13 (0.96-4.69) 


0.031 


2.09 (0.95-4.60) 


0.041 


Cytogenetics 






0.056 








t(8;2l) (24) vs Normal (16) 


0.39 (0. 1 6-0.95) 


0.012 


0.39 (0. 1 6-0.96) 


0.022 




Others (14) vs Normal (16) 


0.54 (0.22-1.32) 


0.165 


0.54 (0.22-1.32) 


0.182 


RI regimen 


IDA + BH-AC (35) vs IDA + Ara-C (14) 


1.32 (0.59-3.01) 


0.512 


1 .47 (0.66-3.29) 


0.373 


BM blast, 7o 


^50 (26) vs <50 (28) 


0.80 (0.39-1.66) 


0.540 


0.71 (0.34-1.48) 


0.347 


Abbreviations: BM 


= bone marrow; CI = confidence inter/al; H = high 


expression; HR = hazard ratio; HSCT = haematopoietic stem cell transplantation; L = low expression; 



RI = remission induction. 
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Figure I Distribution of nnRNA expression levels for 1 2 genes in 54 patients. (A) Schematic flow of therapeutic actions of nucleosides and anthracycline. 
Effector molecules encoded by the genes analysed in this study are marked in bold. (B) The expression of 12 candidate genes evaluated by quantitative 
RT-PCR. Details regarding the processes of raw data are described in the 'Materials and Methods' section. The final data values were expressed on 
loglO-scale x-axis. The solid line in box plots (ranging from 25th to 75th percentile) and the circular bullets indicate the median and outlier values, 
respectively. 
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Figure 2 Kaplan-Meier evaluation of the prognostic impact of RRI\A2, T0P2B and T0P2A for clinical outcomes. According to the expression levels of 
RRM2, T0P2B> and T0P2A, a total of 54 patients were divided into low and high groups. The disease-free sun/ival (A) and overall survival (B) curves were 
generated by Kaplan-Meier survival analyses. Statistical difference (s.d.), HR and P-value were calculated via log-rank test. 



significant correlation of gene markers with cytogenetics was 
observed, and HSCT only showed a weak correlation with T0P2A 
(correlation coefficient, 0.296). 

Next, gene-clustering analysis was conducted in order to 
evaluate the correlations among expression patterns of each 
specific gene. The results visualised a conflict of mRNA expression 



between T0P2B and RRM2 genes. In addition, T0P2B evidenced a 
noticeable trend of counter-expression against its 170kDa iso- 
form-encoding T0P2A gene in the AML samples (Figure 3A). 
Therefore, we attempted to divide the 54 patients into two groups 
according to a reference T0P2B/RRM2 ratio value of 1.0 (range, 
0.015-19.021). Kaplan-Meier estimation curves for DFS and OS 
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Figure 3 Prediction of disease-free survival and overall survival using TOP2B/RRM2 or TOP2B/TOP2A ratios. (A) A 54-patient population (rows) was 
rearranged in accordance with the gene expression profiles of RR/V\2, T0P2B and T0P2A genes (columns) via gene-clustering analysis. (B, C) The patients 
were divided into low and high groups according to the reference value, 1.0, of TOP2B)/RRM2 (B) or TOP2B)/TOP2A (C) ratio. Kaplan-Meier estimation 
created the survival curves for disease-free survival and overall survival. Statistical difference (s.d.), HR and P-value for survival curves were calculated via 
log-rank test. 



showed meaningful differences between both groups (Figure 3B). 
The HRs of the high group against the low group were 0.40 (95% 
CI, 0.19-0.85; P = 0.010) and 0.41 (95% CI, 0.19-0.85; P= 0.014) for 
DPS and OS, respectively, in a univariate analyses conducted using 
log-rank tests. In a similar manner, the total patient sample was 
again dichotomised into two groups by a reference T0P2B/T0P2A 
value of 1.0 (range, 0.005-9.027), allowing both groups to be clearly 
distinguished in Kaplan-Meier Survival estimations (Figure 3C). 
The differences between the two groups subdivided by TOP 2 A 
alone did not reach the threshold for statistical significance 
(Table 2). When the mRNA level for T0P2A was combined with 
that of T0P2B, however, the univariate analysis generated great 
impact and reliable significance for good response with a lowered 
HR of <0.3 (95% CI, 0.14-0.62 for DFS and 0.13-0.60 for OS; 
P< 0.001 for both). 

We further conducted multivariate analyses to explore the 
influence of clinical factors in regard to the prognostic properties 
of four gene markers, which were identified in the univariate 
analyses (Table 3). A Cox proportional hazards model adjusted for 
age, HSCT and cytogenetic abnormalities indicated that a high 
TOP2B/TOP2A ratio was an independent favourable prognostic 
factor for both DFS (HR, 4.17; P = 0.002) and OS (HR, 3.4; 
P = 0.005), as shown in Figure 4, whereas the predictive likelihoods 
of the others fizzled out in this analysis. Cytogenetic status, among 
the clinical features applied, was identified as an independent risk 
factor for DFS (P= 0.020), but not for OS. 

Expression of RRM2y T0P2B and T0P2A in drug- sensitive 
and -resistant AML cell lines 

In an effort to confirm whether the interesting expressions in the 
other candidates as well as the marker genes in the BM blasts 
of patients was also observed in drug-resistant cell models, we 
conducted real-time PCR analysis in Ara-C and idarubicin-resistant 



Table 3 Multivariate analysis of clinical outcome by Cox proportional 
hazards model 



Disease-free survival 



Overall survival 



Variables 



HR (95% CI) P-value HR (95% CI) P-value 



TOP2mOP2A ratio 
High (^1.0, N = 2I) 
Low (< 1.0, N = 33) 

Cytogenetics 

Normal (N= 16) 
t(8;2l) (N = 24) 
Others (N= 14) 



I I 

4. 1 7 ( 1 .73- 1 0. 1 0) 0.002 3.4 ( 1 .5-8. 1 ) 0.005 



0.45 (0.1 7- 1. 1 5) 
1.35 (0.51-3.6!) 



0.035 



0.094 
0.544 



Abbreviations: CI = confidence interval; HR = hazard ratio. Gene nnarkers and 
clinical factors, for example, RRM2, TOP2B, TOP2B/RRM2, TOP2B/TOP2A, age, 
cytogenetics and HSCT, which displayed P-value <0.2 in univariate analysis, were 
considered in this analysis. Exceptively, we included TOP2A in the variables even 
though a predictive significance of this gene was not satisfied to our criteria in 
univariate analysis. Shown are the variables remaining only in a final model of the 
Cox proportional hazard model with a stepwise forward method. 

AML-2 cell lines. The drug-resistant AML-2/IDAC cells evidenced 
higher transcriptional levels of RRM2 and T0P2A than were 
observed in parental AML-2 cells, whereas T0P2B was down- 
regulated in AML-2/IDAC cells. These are quite consistent with 
the results observed in AML patients (Figure 5A). Additionally, 
an unusual downexpression of dCK was showed in AML-2/IDAC 
cells, reconfirming our previous finding that the extremely 
decreased expression of dCK gene is a main reason of Ara-C 
resistance in this cell model (Song et al, 2009b). When verified by 
semi-quantitative PCR and gel electrophoresis, the expression 
patterns of RRM2y T0P2B and T0P2A genes in the cell lines 
were in an agreement with those from real-time PCR (Figure 5B). 
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Figure 4 Survival curves modified by the Cox proportional hazards model. 
A total of 54 patients were divided into low and high groups according to a 
reference TOPIB/TOPIA ratio of 1 .0. The curves for disease-free survival (A) 
and overall survival (B) were rearranged from a multivariate analysis by Cox 
proportional hazards model via the stepwise method. 
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Next, we evaluated the DNA methylation status of the genes for 
RRM2, T0P2B and T0P2A in order to determine whether the 
differently observed mRNA expressions were modulated by some 
epigenetic modification inherent to drug-resistant AML cells. As 
anticipated, we observed that CpG island methylation for RRM2 
and T0P2A was slightly reduced in the AML-2/IDAC cells, but CpG 
island methylation for T0P2B was higher, relative to that observed 
in the wild- type cells (Figure 5C). These results were quite similar 
to the observations noted with the clinical specimens, specifically 
that most patients with poor prognoses evidenced higher RRM2 
and T0P2A, but lower T0P2B. 



DISCUSSION 

Discovering a novel marker that can be applied to the categorisa- 
tion between patient populations with diverse outcomes against 
routine induction regimens at diagnosis is very important. 
Additionally, a more specialised study employing a patient subset 
with a specific subtype is clearly necessary, owing to the high 
heterogeneity of AML. In this study, we searched key determinants 
allowing for the prediction of treatment outcomes in AML patients. 
The results from univariate analysis showed that T0P2B/RRM2 
and T0P2B/T0P2A ratios, as well as individual RRM2 and T0P2B 
expressions, are statistically reliable predictors for DFS and OS. 
When applied to multivariate analysis including other clinical 
factors, as described in Table 3, we concluded that a segregation of 
patients based on the ratios of the two T0P2 isoforms was the 
optimal means of predicting treatment outcomes; this also implies 
an independence of the marker as a prognostic determinant 
against chemotherapy. 

The most interesting finding in this study is that it may prove a 
possibility to predict treatment responses by evaluating T0P2B 
and T0P2A expression levels, showing the conflicting expression 
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Figure 5 Expression and methylation analyses for the marker genes in drug-sensitive and -resistant AML cell lines. (A) mRNA extracted from AML-2 
and AML-2/IDAC cells were applied to semi-quantitative RT-PCR for 12 genes. The results were obtained from the two independent experiments and 
expressed as relative values against mRNA levels in AML-2 cells. N.D, not determined. (B) mRNA extracted from AML-2 (WT) and AML-2/IDAC (IDAC) 
cells were also used for semi-quantitative RT-PCR for RRM2, TOPIB and T0P2A. GAPD\-\ was used as an internal control. Using densitometric analysis, the 
results were expressed as the means ± s.d. from three independent experiments. *P<0.05. (C) gDNA from WT and IDAC cells were applied to MSP 
analyses for unmethylated (UM) or methylated (M) sequences of RRM2, TOPIB and T0P2A. The PCR products were analysed via agarose gel 
electrophoresis. Densitometric analysis was conducted, and the results were expressed as mean ratios ± s.d. between M and UM bands from at least two 
independent experiments. *P<0.05. 
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trends between the two genes. Although we anticipated finding 
overexpression of both genes in the patient group with the shorter 
survival rate, the majority of cases (43 of 54 patients, 79.6%) 
evidenced a differential pattern, with more than a 1.5-fold 
difference in mRNA expression levels noted between the two 
genes. Thus, we concluded that the prediction of DFS rate via the 
ratio of T0P2B/T0P2A as opposed to that of T0P2B alone is more 
accurate, and maximises the difference between the two patient 
populations (12.9 vs 5.2) in the univariate analysis. The T0P2 
isoforms have been identified as essential nuclear enzymes, which 
potentially encourage the rapid replication of tumour cells by 
modifying DNA structures (Vos et a/, 2011), but their expression 
and roles in sensitivity against anthracycline- containing regimens 
are controversial as ever. Most studies have focused on T0P2A and 
reported involvement of tumour chemoresistance to DNA- dama- 
ging agents in AML as well as other tumours (Chekerov et al, 2006; 
Chen et a/, 2009; Coss et a/, 2009). Some groups suggested no 
significant correlation between topoisomerase levels and treatment 
outcome at diagnosis (Kaufmann et al, 1994; McKenna et al, 1994). 
Another recent study indicated that a break of double-strand DNA 
by T0P2B, as opposed to that by T0P2A, is more closely 
associated with the development of secondary malignancy after 
treatment with T0P2-targeting drugs such as etoposide and 
doxorubicin (Azarova et al, 2007). With regard to the case of 
cancer cells that escape from or resist anticancer drugs, one 
possibility is that reduced topoisomerase II levels occasionally 
allow the cells to survive from cytotoxic attack by maintaining 
their quiescence (Friche et al, 1991; Chen et al, 2011). Apart from 
relationship studies between topoisomerase expression and drug 
response, it has been reported that T0P2A activity but not T0P2B 
in AML blast was decreased by treating with daunorubicin and 
idarubicin, indicating that a primary target of anthracycline is 
T0P2A (Gieseler et al, 1996). Microarray-based approach of 
diagnostic AML with normal cytogenetics showed that T0P2B 
expression was much higher in patients showing continuous CR 
than in refractory or relapsed patients after CR, although the 
researchers did not evidence a direct correlation between T0P2B 
expression and sensitivity to the induction chemotherapy (Vey 
et aly 2004). Therefore, our results imply the existence of a 
mechanism whereby high T0P2B assists rapid proliferation of 
leukaemia cells and induces a high degree of susceptibility to 
chemotherapy, thereby indicating that T0P2B as well as T0P2A 
manage in drug sensitivity against anthracyclines by distinct 
independent ways. Another possible explanation is that T0P2B/ 
T0P2A ratio might be not relevant to a susceptibility against 
antitumour agents but to a biological feature of AML with good or 
poor outcome, empathising the value of the parameter as a 
predictive marker only. However, the functional and differential 
roles of both topoisomerases remain to be clearly elucidated. The 
reason for the differential regulation of their expression in AML 
cells will require elucidation in future studies. 

Previous in vivo and in vitro studies have demonstrated that 
nucleoside drug metabolism-related genes may have an important 
role in responses against the drugs and can be used as a prognostic 
indicator not only in leukaemias. One study using samples of 
Ara-C- treated AML patients indicated that high dCK and low 
NT5C expression were noted in the group with longer survival 
rates (Galmarini et al, 2003). Additionally, the roles of dCK, CDA 
and hENTl were addressed in several studies employing in vitro 
leukaemia cell models (Ge et al, 2004; Cai et al, 2008; Song et al, 
2009b). In our study, despite its identification as a pivotal key 
enzyme, dCK was not statistically correlated with DFS (P = 0.965) 
or OS (P = 0.852), according to univariate analyses. We also 
detected no meaningful significance in survival analyses using 
other candidates. The discrepancy between our results and those of 
previous studies may be attributable to the specialised role of dCK 
in the activity of anticancer drugs, and to the fact that other studies 
have employed differently subtyped AML populations. Indeed, to 



screen outcome predictors, other research groups have employed 
mixed populations of AML patients treated with non-unified 
regimens, whereas in this study we used more focused populations 
in the aspect of disease type and chemotherapy regimen. Additionally, 
our previous in vitro study showing that dCK performed a critical 
function in Ara-C resistance, but not in idarubicin resistance, of 
AML cells bolsters this discrepant observation (Song et al, 2009b). 

We also attempted to address the regulation mechanisms 
underlying the differential expression of marker genes between 
AML patients with good or poor outcomes. DNA methylation is a 
well-established mechanism of gene regulation in AML (Figueroa 
et a/, 2010). One study using 42 AML samples suggested that a 
poor prognosis by overexpression of ABCBl was correlated with 
hypomethylation in its promoter region (Nakayama et al, 1998). In 
order to confirm any correlation between methylation status of the 
selected genes and drug resistance, we employed a drug-resistant 
in vitro leukaemia model. As described in a previous report, an 
AML-2/IDAC (double-resistant to Ara-C and idarubicin) cell line 
was originated by way of a chronic drug exposure of the drug- 
sensitive parental AML-2 cells (Song et al, 2009a). We hypothesised 
that RRM2 and T0P2A expression would be downregulated in drug- 
resistant cells, whereas T0P2B expression would be upregulated. 
As anticipated, altered gene expressions were noted in AML-2/IDAC 
cells, as compared with the AML-2 cells (Figures 5 A and B). We also 
determined that the differential expression of the genes can be 
regulated by the DNA methylation status in the CpG islands 
(Figure 5C). Although AML-2 cells are not fully representative of 
the clinical AML blasts used in this study because both resistant 
cell lines were already exposed to the drugs, these results are 
reflective of a connection between the RRM2, T0P2A and T0P2B 
genes with the drug-resistant phenotype in AML cells. Additionally, 
this cell line may prove useful as a good model for in vitro attempts 
to define the contribution of the markers in response to treatment. 

In summary, the results of this study demonstrated that the 
expression ratio of the T0P2B and T0P2A genes was correlated 
with treatment outcomes in AML patients treated with cytosine 
arabinoside-based chemotherapy in an independent manner. The 
finding may prove useful in the prediction of treatment responses 
in M2-subtype AML patients, and also expands the range of possible 
treatment choices for improved therapeutic efficacy. However, addi- 
tional confirmatory studies in validation sets wiU be required in order 
to actualise the possible benefits illustrated by the results of this study. 
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